A model for ensemble NMR quantum computer using 

antiferromagnetic structure 

A. A. Kokin 



Institute of Physics and Technology, Russian Academy of Sciences, Nakhimovskii pr. 34, 
Moscow 117218, Russia, E-mail: kokin@poboxes.com 

Abstract 

The one-dimensional homonuclear periodic array of nuclear spins I = 1/2, owing 
to hyperfine interaction of nuclear spins with electronic magnetic moments in anti- 
ferromagnetic structure, is considered. The neighbor nuclear spins in such array are 
opposite oriented and have resonant frequencies determined by hyperfine interaction 
constant, applied magnetic field value and interaction with the left and right nuclear 
neighbor spins. The resonant frequencies difference of nuclear spins, when the neighbor 
spins have different and the same states, is used to control the spin dynamics by means 
of selective resonant RF-pulses both for single nuclear spins and for ensemble of nuclear 
spins with the same resonant frequency. 

A model for the NMR quantum computer of cellular-automata type based on an 
one-dimensional homonuclear periodic array of spins is proposed. This model may 
be generalized to a large ensemble of parallel working one-dimensional arrays and to 
two-dimensional and three-dimensional structures. 

Introduction 

The fundamental obstacle, preventing experimentalists from extending the number of qubits 
to L > 1 in an individual molecule of the liquid-state NMR quantum computer, is the 
difficulty of distinguishing L unique set of two-state cells. To remove this obstacle it was 
already proposed several models for solid-state quantum computers with both individual 
and ensemble control of qubits. One of such potentially realizable model based on a one- 
dimensional cellular automaton, using an one- dimensional periodic array ABCABC. . . of 
three types of two-state quantum-mechanical cells (they may be heteronuclear system of spins 
/ = 1/2) with distinct resonant frequencies and local interaction between near neighbors, 
was first considered by S.Lloyd |L|. The effect of the interaction contains in a shift of the each 
cell energy levels depending on states of its neighbors. After using the resonant tt— pulse all 
cells of type A, for instance, invert their state if, and only if, the left neighbor C is in ground 
state and the B on its right is in excited state. In M, it was represented algorithm, which 
was applied globally to all cells, so that there is no need to address cells individually. This 
model was recently developed by S.Lloyd f2j. 
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The more general model of a solid state ensemble NMR quantum computer was described 
in [[|, where it was considered periodic structure of ABCABCABC. . . type in two or three 
dimensions with the nuclear spins 1/2 only of three distinguish types A, B, C. It was supposed 
that the nuclei are embedded in a crystal lattice of some solid state compound with spinless 
nuclei and all spins are initialized to the ground state |0). Each ABC-unit of this superlattice 
can be used to store quantum information by setting one of spin up or down. This information 
can be moved around via some quantum cellular state shifting mechanism. Cascading unitary 
quantum SWAP operations of A <^ B, B C, C -v^ A, A <^ B,. . . is used for this process. 
An ancillary dopand nucleus D with spin 1/2 in the proximity of an A-site can serve as 
the input /output port. A local environment region near dopand nucleus provides a large 
quantum system with a wealth of qubits and only three types of nuclear spins. Therefore 
impurity doping may induce large-scale quantum automata in a single crystal and the whole 
crystal contains a huge ensemble of such identical NMR quantum computers — large artificial 
" molecules" . 

One-dimensional scheme that was based only on two different A and B types of cell in a 
periodic array without the ability to distinguish the left neighbor from the right was described 
in ||. Each two-state cell of the scheme has ground | j) and excited | f ) internal eigenstates 
and can represent any quantum superposition of these states. All cells are initially in the 
same ground states | j) and the state of the all array is | HI ... j), similar to an one- 

ABA 

dimensional two-sublattice ferromagnetic. Each qubit of quantum information in the state 
is represented by four consecutive units: the qubit basis state "0" is represented by unit 
I ttll)> whilst the state "1" is represented by | ||Tt)- 

ABAB ABAB 

The model of array described below could be realized by using a linear artificial "molecule" 
with A and B cells alternating along its length in antiferromagnetic-type structure. As the 
cells in this array are used only identical nuclear spins I = 1/2. The neighbor nuclear spins 
in the ground state of antiferromagnetic structure are opposite orientated and have distinct 
resonant frequencies determined by hyperfine interaction constant, by applied magnetic field 
value and by interaction with the left and right nuclear neighbor spins. The major advantage 
of this variant over the ferromagnetic structure is that the antiferromagnet doesn't have the 
total spontaneous magnetization and the nuclear resonance frequency doesn't depend on the 
sample shape. 

1 The one-dimensional antiferromagnetic model on 
atoms 31 P. 

In |5], H it was suggested a bulk-ensemble generalization of the silicon quantum computer 
model proposed by Kane previously @. In ensemble case, unlike the individual Kane's 
model, two-type electrodes A and J form a set of narrow (/a ~ 10 nm) and long (several 
micrometers) strips. The distance between neighbors A gates was assumed Z x ~ I a- Along 
the gates A, donor 31 P atoms l y distant from each other are placed. If exchange interaction 
constant for localized electronic spins along the strip gates is more than for electronic spins 
between neighboring strips and more than Zeeman energy J(/ y ) J(/ x ),2/i B -B (£> is the 
induction of the applied magnetic field), it produces an artificial one- dimensional antiferro- 



2 



magnetically ordered state of electronic spins. At the temperatures well below the critical 
temperature (Neel temperature) Tns ~ J(l y )/k (k — the Boltzmann constant) we will have 
a pure macroscopic electronic ground quantum state. Due to hyperfme interaction nuclear 
spins will be oriented according to the electronic spin direction in the resultant field and 
will form array with the alternating orientation of nuclear spins. Note, this state is not the 
true pure nuclear antiferromagnetic state so as long as the phases of distinct nuclear spins 
at macroscopic distances are not correlated at temperature of order or higher then critical 
temperature of nuclear magnetic dipole ordering, that is T > T NI ~ (1CT 6 — 10~ 7 )-K" ||. 
However, the phase correlations of near neighbor nuclear spins of course exist. 

The nuclear resonant frequencies z^b of neighbor nuclear spins are different for each 
of the magnetic quasi-one-dimensional subarrays A and B in the chain and depend on the 
states of neighboring spins. We will take it in the form: 

^a,bO< + m>) « ± A/2 - J n (m< + m > )\/2irh, (1) 

where /zn = 5.05 • 1CT 27 J/T is the nuclear magneton, A is hyperfme interaction constant, 
(for 31 P: <7n = 2.26, A = 7.76 ■ 10~ 26 J), I n — the constant of two neighbor nuclear indirect 
spin-spin interaction, m < and m> are the magnetic quantum numbers for the left and right 
spins. The nonsecular part of interaction is neglected here taking in to account that g^fi-^B, 
A/2 > J n . 

Thus we have the one-dimensional homonuclear periodic array of nuclear spins I = 
1/2, formed in the one-dimensional antiferromagnet at the applied magnetic field, owing to 
hyperfme interaction of nuclear spins with the electronic magnetic moments. At magnetic 
fields B > A/2g^fi^ ~ 3.5 T and at temperatures T ~ 10~ 3 K the nuclear spins 31 P have in 
each subarray almost 100% orientation {2nhv^/kT < 1), that is they are in ground state. 
Note, that the using of dynamic methods, such as optical pumping, makes possible the high 
orientation of nuclear spins also at more large temperatures. 

We will estimate here the exchange interaction constant J > 2^,-qB ~ 6.5 ■ 10~ 23 J, the 
critical temperature T^s ~ J/k ~ 4.5 K and the nuclear spin critical temperature, that is due 
mainly to the Suhl-Nakamura indirect spin-spin interaction, Tni ~ I n /k ~ A 2 / Jk ~ 10~ 5 K. 

Here we shall use for the organization of logic operations the addressing to spin states and 
qubits, analogously to consideration Q . We shall consider at first the simple one-dimensional 
model of the antiferromagnet, in which each cell is represented by the magnetic atom and has 
one electronic and one nuclear spin 1/2 with hyperfme interaction, similar to the mentioned 
above artificial molecule of antiferromagnetically ordered donors 31 P in silicon substrate. 

Nuclear spins of identical atoms at g^fi^B < A/2 are oriented according to the electronic 
spin direction in the resultant field and will form a periodic ground state array of ABAB. . . 
type: TIT I • • •> where j marks the ground state of nuclear spin in an A-site and J. — the 
ground state of nuclear spin in a B-site, that is we have here homonuclear system of spins 
at two distinct ground states. Each nuclear spin in A-site of this scheme, has two internal 
eigenstates — ground | T) and excited | JJ.) and in B-site, accordingly, - | |) and | ff). We 
take into account that the life time (the longitudinal or spin-lattice relaxation time 7\) of 
excited states at low temperatures is very long. Each qubit of quantum information in this 
state will be represented here, similar to Q, by the four consecutive cells: the logical qubit 
basis state "0" will be represented by unit | -ILftTI)) whilst the state "1" — by | TI4lT)- It 
is important here that the resonant frequencies of nuclear spins depend on neighbor spins 
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states. 

The input and output of the information in the array of ground states spins could be 
performed at the ends of the array, where the nuclear spin (say in A-site at the left end) has 
only one neighbor spin and distinguishing resonant frequency Vx-x/2 (m< + m > — — 1/2)- 
The corresponding selective resonance RF 7rA,-i/2 — pulse inverts only one nuclear spin (in 
A-site) at the end of array and doesn't influence on any ones. Then the new selective RF 
^3,0" pulse will invert next nuclear spin (in B-site), which has the opposite orientation of 
ground and exited neighbor nuclear spins (m< + m > = in A-site) and consequently the 
new resonant frequency, distinguished from the frequency of spins with the neighbor nuclear 
spin in ground states (m< + m > = 1). Thus the qubit state "0", that is | -lift? I)? is formed 
in the following way (the pulses act on underlined spins): 



"0" 

7r A,-i/2~P u l sc ,, ... . 7T B ,o-pulse : 



nut... ■ ^ writ... wit... 

-BAB A-AB ABAB 

The qubit state " 1" at the edge of array is formed by means of still three pulses: at first 
7Ta,o, then 7Ta,-i/2 and 7r B ,o — pulses: 

tta.,0 — pulse ta,-i/2— pulse 



MTIT... ilrUT.-- ^ ttllt 

AB-B - - 



""B.o-pulse 



ABAB 



The states |fHJ4t) and HT1HI-) m ay be called as the reversed states relative to the states 



\Un) and IMTI). 

Note that a random inversion of only one spin will result in completely destruction of the 
qubit. However, to form, for example, the error of "0" =>- "1" type in the coding of stored 
quantum information it is essential to invert simultaneously four spins. Therefore, it may 
be concluded that the considered way of qubit coding ensures a better fault-tolerance with 
respect to this type of errors. 

Authors of have considered also another scheme of the four-spin encoding two logical 
qubits, which are represented by the two zero-total states of four spins, generated by the pairs 
respectively of the singlet and triplet states. This scheme leads in the collective decoherence 
conditions to the fault-tolerant implementation of quantum computations. The collective 
decoherence conditions can be attained in coupled spins at very low temperatures, where all 
collective but the longest wavelength acoustic phonon modes are quenched. 

The further shift-loading of qubit states into the array is implemented by means of pulse 
sequence 7Ta,o, ^b,o, tta.o, ttb.o ■ - •, which is represented by following SWAP operation: 



ft A 0~ pulse 7Tb q— pulse 7ta,o — pulse 



ab-b 

i"B,o-pulse 



tit i m I--- tmuMuu- 

- - AB AB 



and so on. 
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The role of the atoms at array ends can play here, as it was discussed in ||, also dopand 
nuclei D at the certain place of the array with distinct resonant frequency or a defect that 
modifies the resonant frequency of the nearest nuclear spin in the array. Starting from the 
perfectly initialized states inputting the information can be performed by setting the dopant 
D-spin to desired state by means of pulse at his resonant frequency. The nuclear spin state 
of cell nearest to the dopand is created by SWAP operation mentioned above. After the 
required information is loaded, D-spin is reset to the ground state |0). Upon completion of 
computation, the state of any qubits can be measured by moving it to the A-site nearest to 
D, then swapping A <^ D and finally measuring the state of D-spin. 

2 One-qubit operations in one dimension 

As in H], we introduce still 7Ta,i and 7Tb,-i— pulses and operators Ua,i and Ub,-i- The last 
means, that each spin in A- and B-site is subjected to a unitary transform U, which acts on 
spins in A- and B-sites with resonant frequency corresponding to the both neighbor nuclear 
spin in the same excited states (m < + m > = ±1, sign "+" is for excited neighbor spins in 
B-sites | ft), and sign "— " — for A-sites | JJ-), see Table). Note, that when operator U is a 
simple inversion, the actions of Ua,i and Ub,-i are identical to 7Ta,i and 7^,-1— pulse. 

Table. The tt— pulses for spins in A- and B-sites 



Neighbor spin states. 


A-site 


u 


m 


in 


im 


tn 














A 


Resonance frequency 




M-l/2) 


"a(-1) 


^a(O) 


^a(O) 


"a(1) 


n— pulses 




7TA-1/2 


tta -1 


7TA,0 


tta.o 


7TA,1 


Neighbor spin states. 


B-site 


IT 


TIT 




U1 


m 








B 






B 


Resonance frequency 




u B (l/2) 


Mi) 




Mo) 


^b(-I) 


tt— pulses 




7Tb,1/2 


7Tb,1 


7I"B,0 


7TB,0 


7TB.-1 



Let us construct now the logical gates of quantum computer. At first we shall investigate 
at first the scheme for one-qubit gate. The considered section of the array (Fig. 1) contains 
three qubits in states "1", "0" and "1", each being separated by number multiple four of 
spacer cells -nuclear spins in ground state. Therefore each qubit requires a total of eight 
physical spins in the array (four for the encoding plus four spacers) (Fig. 1). As it was made 
in [1], U , we will also introduce here one control unit, (not to be confused with control qubit 
in case of CNOT gates), which is represented here by six consecutive cells in the pattern 
fNJ-ITfNJ-- The control unit (CU) exists only in one place along the array and is separated by 
odd number of spacer cells — spins (the scheme at Fig. 1 has three). 
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"0" 


"1" 






cu 






i 1 1 1 1 




1 1 

ii i 


1 1 1 




T =* 


n =*. . . t 

A, 1 


If til 


t i 1 1 1 1 1 


t j 


lit 




J =* 


n =>. . . t 

B, 1 


lUjt 


t i t i 1 1 1 


1 1 


1 J 1 


I • i , . , 


J =* 


n =>. . . II 

A, f 


1 1 1 1 1 




u 


1 J 1 
cu 


nut j 


J =* 


n =*. . . II 

B, t 


1 1 i t i 


1 1 1 1 1 1 t 


tj 


1 1 1 


IIUU 


J =* 


n =*. . . t 

A, o 1 


t II 1 1 


1 1 u t n 


1 1 


ill 


1 j U tJ 


J =* 


B, 1 


1 1 1 1 i 




it 


t Ii 


ptltl 


J =* 


n =*. . . T 

A, 1 


it in 
■■I" 


H H H 


u 


1 1 1 


turn 

"1" 


J =* 


n =>. . . T 

E, 1 




t i 1 1 1 1 1 


tj 


i t i 




T ' ^ 


n =>. . . t 

A, o 1 


1 1 1 1 1 


u 1 1 1 1 1 

K K K K 


1 1 


1 1 1 


Uttii 

"1" 


T ' ^ 


B, 1 


1 U U 


until 


1 1 


1 1 1 




j =* 


A, 1 


J U U 


mutt 


tl 


i 1 1 




1 



vo-'-U 1 1 n u u u t n 1 1 j t i t j t n* 

v. *■ ■ ■ t i t i 07117 1 Hit jTTT 1 1 1 1 n * 
v.*---t i n n 1 1 u "ft 1 1 1 1 jTT] 1 1 tit* 



v.'-t n 1 1 n 1 1 1 1 1 n 1 1 1 1 n 1 1 1 1 1* 

CU "1" "O" 

v.*-t t ii 1 1 u 1 1 1 n 1 1 1 1 i 1 1 1 i t i t ■ 



Fig. 1. The scheme of the SWAP update pulse sequence. 



The applying SWAP update sequence of pulses 7Ta,o, ^b,o, tta.q, ^b,o, ^a,o ■ ■ ■ moves the 
qubits to the right and CU to the left relative to the qubits, yet the form of the qubits 
and the CU are preserved. The CU passes through qubit in state "1" and "0", leaving it 
unchanged and continues further (Fig. 1). 

To implement the one-qubit logical gate the additional computing updates six pulse 
sequence 7Tai, ttbi, ttbo, ^ai, ^bo, Ua,i is applied when CU reaches the mid-way through 
passing the qubit Y ("1" or "0", marked * and ** at Fig. 1). The effect of the additional 
sequence is to apply a unitary transform U only to the spin representing the qubit Y: t = UY 
(Fig. 2). The scheme of the additional computing update pulse sequence after stage * is 
shown at Fig. 2. The scheme of sequence after stage ** is shown in Appendix Al. 



v.*-l 1 t , 


i 1 1 1 1 1 1 1 1 


t n n n 1 1 1 1 1 1 


Vi'-Ht, 


i 1 1 Ii t i 1 1 


tnitntnnt 


•..i*-»lt. 


ittut Jtl 

"O" 


t n 1 1 1 \ 1 1 i t i t 


V. 4 -! 1 1 j 


, t n 1 1 1 1 1 


1 1 1 ! J (1 1 i 1 1 1 1 1 


1 1 J 


, t n 1 1 1 1 1 


1 1 1 1 1 m 1 1 1 1 1 1 
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vo*---i 1 1 i r t n 1 1 1 1 1 1 1 1 1 n 1 1 i t i t* 
w A .i *-i 1 1 1 1 TTu 1 1 1 1 m 1 1 n't 1 1 1 1 1 . 

Fig. 2. The scheme of the additional computing update pulse sequence after stage * 
Let us take now the one-cell unitary operation at the end of additional sequence: 

U A ,il T> = a| T) + &U),H 2 + |&| 2 = 1. (2) 

Re-applying the update pulses after unitary operation in reverse order the CU moves 
away from transformed cell and is returned to its initial state (see Appendix A2). We will 
have in the result the superposition of states: 

"1" "0" 

\^) = a\ nU) + b\UTi)- (3) 

The CU and additional computing updates pulse sequence together ensure the computing 
operations with qubits. Note, such one-qubit gate requires seventeen elementary pulses. 

Not let us consider the state of quantum register, shown at the first line on Fig. 1 and 
apply after stage marked ** the pulse 7Ta,i. The result is presented on Fig. 3. 



w-UUtttnttttttitiutti jn 



Vi'-t 1 1 | 1 1 n | t | t | t t | t | 1 | 1 1 ( | t* 

Fig. 3. The result of the SWAP pulse sequence finished by pulse 7Ta,i- 

We see that the CU moves transparently past the qubit "1" and continues until mid- way 
through passing qubit "0". Now the CU itself is subject to a transformation: it is altered 

from lULitlUL to llTllTIlT (o m y for passing the qubit "0"!) and qubit "0" itself will be 
destroyed (Fig. 3). 

Now we will apply again the SWAP pulse sequence and will become: 
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b, o 



a, o 



B, 



A, 



B, 



n . . 

A, 



B, 



a, o 



1 1 ninn 



1 1(| 1 1 ( 1 1 

t tTtTjl 1 1 
jTftlTi 1 t 




( t ( t t | t t | 



t j t n t i t n 



1 1 1 1 1 1 1 1 1 1 



Fig. 4. The scheme of the SWAP pulse sequence following sequence Fig. 3. 
We see here, that by passing qubit "1" the altered CU is preserved his form. 

3 Two-qubit operations in one dimension 

To implement the two-qubit gate such as CNOT it should be applied the another update 
pulse sequence. 

Let now the qubit "0" will be as control qubit of CNOT gate. The CU transforms in 
altered form by passing the control qubit and then we extend the sequence Fig. 5 after stage 
marked ***** by following pulses with the end inversion pulse Ua,i = tta,i : 



KA,1, 7TB.1, VT B ,o, 7T A>0 , 7T A ,l, VT Bi0 , 7T A ,o, 7T B ,l, ^A,l- 



(4) 



We will obtain: 
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A, 



n . . 

A j 1 



B, 1 



B, 



a, o 



A, 1 



B, 



A, 



B, 1 



A, 1 



ft 



t Ji J 



tin 
t j t j 



nil 



i u t 
i . i ' 



1 1 1 1 



t j t j 



nil 



nil 



Fig. 5. The scheme of the update pulse sequence (H) 

The last inversion operation doesn't have effect on target qubit "1", as it must be for 
CNOT gate. The reverse sequence returns CU and qubits to their initial states. 

Let us return then to sequence Fig. 1 and continue it after stages marked **, when CU 
passes mid-way through qubit "1" by following sequence: 



a, o 



1 1 Jl 1 1 1 1 U 1 1J 1 1 1J n 1 1 1 1 1 1 



an 



Vi-*-T t J 1 r 1 1 1 u 1 1 J 1 1 1 J 1 U 1 1 1 1 N 



Vi--t t J 1 1 1 1 1 u 1 1 J 1 1 1 j t u 1 1 1 1 1* 
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B, 



A, 



A j 1 



B, 



a, o 



B, 1 



A, 1 





Fig. 6. The scheme of the update pulse sequences after stage ** when CU passes the 
qubit " 1" 



The scheme of the update pulse sequences after stage **** when CU also passes the qubit 
" 1" is shown in Appendix A3. 

We see, that the last inversion operation has effect on target qubit, as it must be for 
CNOT gate only when the control qubit is " 1" . 

We can consider a large set of quasi-one-dimensional antiferromagnetically ordered weakly 
coupled at J(l y ) 3> J{l x ) arrays donors 31 P in silicon substrate as an ensemble of artificial 
molecules. In this case, there is no need to address qubits individually. We suppose that for 
the determination of nuclear spin states in ensemble of those identical artificial molecules, 
as in liquid-state bulk-ensemble quantum computer ||, [10| , there are no need to fulfill the 
electrical measurements and consequently any electrodes. Since the read-out signal in this 
case will be proportional to the numbers of artificial molecule in the ensemble, it may be 
used the NMR or fluorescence techniques for ensemble measurement of spin states. 



4 Two- and three-dimensional antiferromagnetic struc- 
tures 

Instead of generalization to the parallel model employing "sub-computers" with one- dimensional 
structure, which was considered in ||, our approach allows to use also two and three- 
dimensional structures. The coupled antiferromagnetically ordered chains model can be 
extended to a two-dimensional antiferromagnetic chess-type ordering. Let the electronic 
spins of the neighboring chains are setting for J(l y ) ^ J(l x ) > 2fi B B. The electronic spins of 
two neighbor chains will be in the singlet ground state. The subarray of nuclear spins will 
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have the opposite orientation of nuclear spins relative to the subarray of neighboring chain. 
The electron subsystem of two neighbor chains is in antiphase state, that is have the half 
period shift of antiferromagnetically ordered electronic spins in the one chain relative to the 
other. The nuclear subsystem of the both chains becomes corresponding chess-type ordering 
(Fig. 7). 

Let us suppose that an initial state containing some number of cells is loaded in the two- 
dimensional structure of many coupled chains. The inputting the information into the cell 
nearest to dopant atoms D (D-spin) is performed by means of corresponding 7r— pulse: |=Hf 
or T=^-U- Then, the resonant frequencies of neighbor spins in the same chain and of spins 
in neighbor chain is altered and another it— pulse will invert one of they or both according 
to the values I\ and I2 of indirect nuclear spins interaction inside and between the chains. 
Therefore, the excited nuclear spin states and accordingly the qubits may be passed to any 
place in all two-dimensional structure. We will suppose as before that the qubits state will be 
represented by four spin states in chains. The computing operation can be fulfilled analogous 
to the above- considered one-dimensional scheme. 



D I "1" 

© 1 _ 



«t«|»f»|»t»J»t J J |»|»|»f»J» 



V I 

Z 



«J«f»|»t»|»t J J |»f» l»f» 

V V V V 

V V V V 

V V "0" V V 



«t«l»t»l»l t t J»t»J»t»J»t J I t 



Fig. 7. The scheme of two-dimensional nuclear spin ordering in antiferromagnetic struc- 
ture. It is showed the different ways that connect the D-spin (marked 0) and a certain 
qubit. 



For everyone CU we have a two-dimensional section of array with large enough number 
of qubits and one dopant atom. It is defined as a single-domain antiferromagnetic sample. 
There are many ways that connect the D-spin and the qubit (Fig. 7). This section plays role 
of many-qubit artificial molecule and the whole structure represent a large ensemble of such 
molecules, which work simultaneously and ensure the parallelism of quantum operations. 

The structures with two and three-dimensional antiferromagnetic and ferrimagnetic or- 
der may be found perhaps among the natural rare earth or transition element dielectric 
compounds. 

The electron magnetization of a one subarray of antiferromagnet is defined by expression 
0: 

2(i B N(S iz ) = 2/i B JV(l - P{T) - V)/2, (5) 
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where for low temperatures in spin-wave approximation 



P ^ (2vr) d / exp(e(A;) jkT - 1 ^ 



is the contribution of thermal and ip — of quantum fluctuations, e(k) = ye^ + (Jak) 2 - 
spin-wave spectrum, a — lattice period, Z — number of near neighbors, d — dimension of 
structure. 

For antiferromagnetic state of easy-axis type, when the interaction Hamiltonian of two 
electron spin j and g for single-axis crystal has the form: 

where J > Ja > 0, Ja — anisotropy constant, the contribution of quantum fluctuation, as 
shown in ||12| , ip = 0. In addition, for k =>- 0, T <C e , where eo ~ ZV J J a, 



P{T) ~ Const • {kTe /J 2 ) d/2 exp{-e /kT) 0, (8) 

that is in this state the thermal fluctuations in electron system also are not essential for the 
ground electron spin states. 

Note, that in the case of easy-flat state, ip ^ 0, but the NMR resonance frequency depends 
on the neighbor nuclear states to a greater extent than in case of easy-axis state |13| . 

The quantum state decoherence at low temperatures is defined on the one hand by the 
active role of electron spin- wave effects They generate the fluctuated local field due 

to Raman process of the electron spin wave scattering on individual nuclear spins. The 
decoherence time or transverse relaxation time T2 of NMR in antiferromagnet for the low 
temperatures (e /kT ^> 1) then is determined by expression ||13| 



1/T 2 ~ {A 2 / J) ■ (kT/J) 3 (e /kT) exp(-e /kT)/n 2 h 0, (9) 

value T 2 rapidly grows. 

On the other hand decoherence is defined by inhomogeneity of the local magnetic fields 
and spread in resonance frequencies. The nuclear spin-spin interaction in natural dielectric 
antiferromagnets is defined mainly by the Suhl-Nakamura indirect mechanism of interaction 
through exchange of spin waves and is typically greater than the value, determined by the 
direct nuclear spin-spin dipole interaction. This interaction of nuclear spins could play a 
large role in the case of high spin concentration. Both of these decoherence mechanisms 
can be, in principle, suppressed by some NMR many-pulse methods using the stroboscopic 
observation of spin dynamics JT4], [L5[] . 

The general requirements for natural antiferromagnetic structures, required for the con- 
struction of NMR quantum computers, can be formulated in the following way: 

1) The operating temperature T must correspond to the fully ordered antiferromagnet 
2ns ^ T ^> T NI and to fully polarized nuclear spins T NS A/J ~ A/k > T ^> T NI . From 
where we will have the value T > 10~ 3 K. 

2) The two-dimensional and tree-dimensional magnetic structure must have chess-type 
order (see Fig. 7). 
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3) The magnetic structure must have the easy-axis state of antiferromagnetism in single- 
axis crystals. 

4) The atoms must have nuclear spins / = 1/2. Electron spins may be S < 1/2. 
There are the rare earth compounds of unique thulium stable isotope 169 Tm, that has 

nuclear spin / = 1/2, = 0.458 and makes up 100% of naturally occurring elements with 
stable spinless isotopes of other elements. They can be: Tm 2 3 , TmSi 2 , TmGe 2 and TmSe 
|T7| |. The ground electronic state of magnetic ions Tm 3+ corresponds to S = 1. The 
natural elements O, Si, Ge and Se have, accordingly, nuclear spin containing isotopes (in 
brackets it is shown the isotope occurrence) 17 I = 5/2 (0.04%), 29 Si / = 1/2 (4.7%), 73 Ge 
/ = 9/2 (7.76%), 77 Se / = 1/2 (7.78%). For Yb 2 3 T NS ~ 2.3 K, isotope 171 Yb (14.31%) 
has 1=1/2 and S' = 1/2 (ground state is Kramers doublet). It is known, that compound 
TmSe has the critical temperature for antiferromagnetic transition Tns ~ 2if |[L7| . The 
antiferromagnets with two different nuclear spin I = 1/2, for example FeF 2 with rutile-type 
and TmAg with CsCl— type structure, which have critical temperature 79 K and 9.5 K, may 
be also of interest to the considered questions. Isotopes 57 Fe (2.19%), 19 F (100%), 107 + 109 Ag 
(100%) have according values #n = 0.182, g^ = 5.26 and #n = 0.24. 

In conclusion, we will point out the several advantages of the considered model: it uses 
the antiferromagnetic structure containing only one type of atoms with nuclear spin 1/2, 
it is not needed to have any gate electrodes, the decoherence associated with noise voltage 
is absent, the considered way of qubit coding ensures a better fault-tolerance with respect 
to the generation of wrong qubits, the model admits an ensemble address qubits, it may be 
used as base for development of bulk-ensembles three-dimensional solid-state NMR quantum 
computer. 

The author is grateful to K.A.Valiev for critical reading of the article and useful remarks 
and V.A.Kokin for the help in preparation of this text. 



Appendixes 

Al. The scheme of the additional computing update pulse sequence after stage 
** at Fig. 1: 



\. --t|ttUtjtttttutUtt|lUN 

HH 

vi *-t jTR 1 1 i t uTffft f |TT| i 1 1 1 * 
Vi *-"t 1 1 1 1 i t i t uTffft 1 j 1 1 1 1 1 1 1 * 
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v.-- 1 I t | J 1 1 i t i UtTj 4 T 4 1 i 1 1 T i t - 



Vi--t 1 1 1 1 1 1 1 1 1 U 1 1 U 1 4 1 i 1 1 1 1 1* 



B, 



t i 1 1 1 i t j 1 1 j 1 1 n t u 1 1 n j t j t 



u ■=►. 

A, 1 



1 1 1 1 1 i t i t H 1 1 j j 1 1 1 1 1 j 1 1 1 f 



A2. The scheme of the reverse update pulse sequence after one-qubit operation 

U Ail : 



U * at 

A, 1 1 < 


1 1 (1 1 (1 1 i 1 1 1 1 


t **t i ■ 


ftit- 


f i V 1 1 1 1 


B, 1 i 


1 1 (1 1 (1 1 1 u 1 1 


t **t i ' 


fin* 


f t n t i t 


A, 1 1 i 


[ t (1 J (1 1 1 1 1 1 1 


t **t r 


fur 


f unit 


B, 1 i 


1 1 i J (1 1 1 1 1 1 i 


t '*t i ■ 


ft it- 


f n 1 1 1 1 




\ t i I i U M t i ' 


t **t i ' 


(tif 


f u 1 1 1 1 




1 1 J I J U M U ' 


t **t i ' 


( 1 1 f 


[nit 1 1 


n =* at 

A, 1 i 


\ tii i «n Mr 


t **t i i 


1 1 1 1 1 


1 1 1 1 u t 


n „ „=* ff t 
B, 1 i 


\ 1 1 ii 1 1 1 i it ' 


t **t i i 


Iff 


1 1 1 j 1 1 1 


m . „=* ff t 
A, 1 i 


i n 1 1 j 1 1 j j t fi i 


i **i i ■ 


fltr 


f 1 1 1 1 1 1 


— J- T i 


" t i 1 1 * * j 1 1 1 


titi 


U 1 1 1 


1 1 1 1 1 1 ■ 
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A3. The scheme of the update pulse sequences when CU passes the qubit "1" 
after stage ****: 



A, 



A, 1 



B, 1 



B, 



A, 



n =*. . . 

A, 1 



B, 



A, 



B, 1 



A, 1 



t|||lt 



I II 1 1 



|t|H 



J t J H 



TttjJ 
t *t|. 



1 1 1 1 J 

fTfft 
jTft 1 

J t J 1 1 



1 1 



It 
(1 

1 1 



I i 
U 

"ft 



1 1 



1 1 



1 1 J 
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